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Aim of the study: To perform a system- 
atic review and formal meta-analysis of 
the literature reporting on HPV detection 
in bronchial squamous cell papillomas 
(SCP). 

Material and methods: The literature was 
searched up to June 2012. The effect size 
was calculated as event rate (95% CI), 
with homogeneity testing using Coch- 
ran's Q and P statistics. Meta-regression 
was used to test the impact of study- 
level covariates (HPV detection method, 
geographic origin) on effect size, and 
potential publication bias was estimat- 
ed using funnel plot symmetry 
Results: Fifteen studies were eligible, 
covering 89 bronchial SCPs from differ- 
ent geographic regions. Altogether, 
38 (42.7%) cases tested HPV-positive; 
effect size 0.422 (95% CI: 0.311-0.542; 
fixed effects model), and 0.495 (95% CI: 
0.316-0.675; random effects model). In 
meta-analysis stratified by i) HPV detec- 
tion technique and ii) geographic study 
origin, the between-study heterogene- 
ity was not significant for either; p = 
= 0.348, and p = 0.792, respectively In 
maximum likelihood meta-regression, 
HPV detection method (p = 0.150) and 
geographic origin of the study (p = 0.164) 
were not significant study-level covari- 
ates. Some evidence for publication bias 
was found only among in situ hybridiza- 
tion (ISH)-based studies and among 
studies from Europe, but with a negligible 
effect on summary effect size esti- 
mates. In sensitivity analysis, the meta- 
analytic results were robust to all one- 
by-one study removals. 
Conclusions: In formal meta-regression, 
the variability in HPV detection rates 
reported in bronchial SCPs is not ex- 
plained by the HPV detection method or 
geographic origin of the study 
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Introduction 

The first reports on solitary bronchial squamous cell papilloma (SCP) appeared 

in 1954, when three authors provided descriptions of this new entity [1-3]. It 
was soon recognized that bronchial SCP can also appear as multiple lesions 
called (tracheo)bronchial papillomatosis [4, 5]. Until the mid 1980s solitary 
bronchial SCP was regarded as an extreme rarity [5, 6]. When first reviewed 
in 1985, Maxwell et at. could collect only 11 such cases in the English litera- 
ture [7]. As emphasized before [8], however, this number is a clear underes- 
timate, because in addition to the cases included in their review [5-7, 9-13], 
a large number of bronchial SCPs had been reported, but mostly in non-Eng- 
lish literature [14-53]. Interestingly, this same trend seems to continue in more 
recent reports [54-75]. According to our best estimates, the total number of 
bronchial SCPs (and papillomatosis) reported by now reaches several hundred 
cases [8]. 

Bronchial SCP must be considered a rare disease, and no estimates are avail- 
able on their incidence rates [8]. There are two peculiar features in the clin- 
ical history of these lesions: i) a tendency to spread to multiple sites within 
the bronchial tree [4, 10, 14, 19, 25, 30, 32, 34, 43, 47, 48, 52, 63, 75], and ii) 
a substantial potential for malignant transformation [2, 8, 11, 27, 39, 40, 50, 
51, 57, 58, 66, 71]. These two characteristics strongly implicate an infectious 
etiology of bronchial SCPs. Indeed, their clinical behavior closely resembles 
that of recurrent respiratory papillomatosis (RRP) [8], and natural history of 
sinonasal papillomas. 

While the etiologic role of human papillomavirus (HPV) in RRP is well estab- 
lished [8], the possible involvement of this virus in sinonasal papillomas was 
first suggested by us in 1983 [76], based on detection of HPV antigen expres- 
sion in a single papilloma. This was soon confirmed by demonstration of HPV 
DNA in a series of cases [77], and HPV involvement in sinonasal papillomas 
is well established at present [78]. 

The first clue on possible HPV involvement in bronchial SCPs was provid- 
ed by two cytological reports from 1975 and 1979, when Roglic et al. [79] and 
Rubel et al. [80], respectively, described characteristic cytopathic changes of 
HPV (i.e., koilocytosis) in the sputum samples, thus emphasizing their simi- 
larity to genital HPV lesions (condyloma). This coincides with the observations 
of us, who in 1979 were the first to describe similar evidence in malignant 
bronchial squamous cell neoplasia [81, 82]. However, it took almost a decade 
until Trillo and Guha (in 1988) first demonstrated HPV antigens and identified 
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HPV particles in two broncliial SCPs [83]. Tine interest in bronchial 
SCP as a potential HPV-associated benign tumor increased 
slowly throughout the 1990s [8]. When first reviewed in 2002, 
the literature reporting HPV detection in bronchial SCPs includ- 
ed 9 studies covering 61 SCP cases, of which 29 (47.5%) test- 
ed HPV-positive [84]. 

With the widespread implementation of prophylactic HPV 
vaccines, it has become increasingly important to estimate 
the total disease burden (in addition to genital lesions) poten- 
tially preventable by the current and new generation HPV vac- 
cines [85, 86]. To cast further light on the controversial role 
of HPV in bronchial SCP, it was felt appropriate to conduct 
a systematic review and formal meta-analysis, covering all 
the published literature without any restrictions to the HPV 
detection methods or geographic origin of the study. 

Material and methods 

Data abstraction 

We identified eligible studies by searching MEDLINE (via 
PubMed) and reference lists from original articles, book chap- 
ters and other reviews until June 2012. No language or date- 
of-publication limitations were imposed. The search terms 
included: papillomavirus, HPV, papilloma, solitary, squamous 
cell, bronchus, bronchial, papillomatosis, malignant trans- 
formation. We considered all publications appearing in 
peer-reviewed journals eligible, irrespective of which method 
(see later) was used for HPV detection, provided that the 
report included exact numbers of analyzed cases and those 
testing HPV-positive, necessary for calculation of the event 
rates (= HPV prevalence) and their 95% confidence intervals 
(95% CI). 

Altogether, > 400 abstracts were derived from the data- 
base, covering the years 1954 to 2012. For the present meta- 
analysis, a total of 15 original studies were determined eli- 
gible, all including cases of bronchial papillomas analyzed for 
HPV detection. Because they were included in a recent meta- 
analysis, all studies reporting only bronchial carcinomas were 
excluded from this meta-analysis. 

From the summaries and/or body texts of each eligible 
study, we abstracted the following information: HPV detec- 
tion method, geographic region of the study HPV genotypes 
analyzed and/or detected, total number of cases analyzed, 
number testing HPV-positive, per cent HPV-positivity, 
authors, and publication year. Only the studies reporting HPV 
in solitary bronchial papillomas or cases of bronchial papil- 
lomatosis were included, whereas all studies describing tra- 
cheobronchial spread of RRP and/or their malignant con- 
version were excluded. 

Statistical analyses 

The specific software Comprehensive Meta Analysis™ (Ver- 
sion 2.2.064; Biostat Inc., Englewood, NJ, USA) was used to 
perform the meta-analysis. The software calculates the event 
rates (logit event rates, SE and variance) based on the events 
and sample size data. To assess overall heterogeneity in the 
event rates between the different studies, Cochran's 0 (two- 
sided) homogeneity p value as well as 1^ statistics (for per- 
centage of variation) were used [87]. To evaluate the possible 



publication bias, funnel plots were drawn by plotting the log- 
it event rates by their precision (1/SE) [88]. Funnel plots were 

evaluated for asymmetry usingthe following statistics: i) Begg 
and Mazumdar rank correlation [89], ii) Egger's test of the 
intercept (regression) [90], and iii) Duval and Tweedie's "trim 
and fill" method [91], which imputes the results that are hypo- 
thetically missing due to the publication bias. 

To assess the variation in the event rates (i.e., HPV preva- 
lence) due to the differences between the individual stud- 
ies, we evaluated the study characteristics using stratified 
random-effects meta-analysis and restricted maximum like- 
lihood meta-regression. Stratified meta-analysis allows 
descriptive comparison of the summary event rates across 
the different categories of specified study characteristics. 
Restricted maximum likelihood meta-regression formally 
compares these differences in event rates across the se- 
lected study-level covariates and estimates the among-study 
variance [92]. Given the inherent differences in analytical 
sensitivities between the different HPV detection methods 
- immunohistochemistry (IHC), in situ hybridization (ISH), 
and polymerase chain reaction (PCR) - meta-analyses 
were performed across these strata. Similarly to distinguish 
true study-specific effects from random variation, all analy- 
ses were also stratified by the geographic regions of their 
origin, blamed as one of the reasons for variability in HPV 
prevalence [8, 45]. HPV detection method and geographic 
study origin were also tested as study-level covariates in for- 
mal meta-regression. 

Sensitivity analysis was performed to assess the influence 
of each individual study on the strength and stability of the 
meta-analytic results. Sensitivity analysis runs the analysis 
k {n = 14) times, each time removing one study to show that 
study's impact on the combined effect size. The sensitivity 
of the results to these one-by-one study removals was eval- 
uated by descriptively comparingthe magnitude and precision 
of the random-effects summary event rates (point estimates). 

Results 

Eligible studies 

A total of 15 studies were considered eligible for the pre- 
sent analysis [58, 70, 74, 83, 93-103], comprising 89 bron- 
chial SCPs analyzed by different HPV detection methods. 
Included are both case reports and larger series, comprising 
up to 31 cases analyzed by PCR [98] and 15 papillomas exam- 
ined by ISH [94] (Table 1). The methods used to evaluate the 
HPV involvement include the following: IHC [83], ISH [58, 70, 
74, 93-97, 102], and PCR [98-101, 103]. Based on the avail- 
able data on geographic regions with different HPV preva- 
lence, the studies were categorized into the following 
regions of ohgin: Other Asia (China excluded), Europe, and 
North America. These 15 studies comprise the target of this 
meta-analysis. Of all analyzed 89 bronchial SCPs, 38 (42.7%) 
tested HPV-positive. 

Analytical results 

Point estimates of event rates 

The crude HPV-positivity (38/89) translates to event 
rates (i.e., effect size, summary HPV prevalence) of 0.422 (95% 
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Table 1. Studies reporting on HPV detection in solitary bronchial squamous cell papillomas 



Method/histological type 






MrV positive 








Detection 


Histological 


j\fe^ _ 




Ml 1 inKpr/ 

INLII 1 lUCI / 


% 


MLlLIIUlo allU yCdi 


Ref. No 


IHC, EM 


SCP 


Canada 


- 


2/2 


100 


Trillo ef al. 1988 


[83] 














Kerley etai 1989 


[93] j 


ISH 


SCP 


France 


11 


1/1 


100 


Bejui-Thivolet et at. 1990 


[58] 




SCP 


UK 


16 


1/15 


6.7 


Carey et at. 1990 


[94] ^ 


ISH 


SCP 


Austria 


6/11 


6/6 


100 


Popper at. 1992 


[95] 


ISH 


SCP 


USA 


6/11 


2/2 


100 


Yousem et at. 1992 


[96] 


ISH 


SCP 


USA 


6/11 


1/1 


100 


Katial etai 1994 


[97] 


r PCR, ISH 


SCP 


Austria 


6, 11, 16, 18 


11/31 


35.5 


Popper et at. 1994 


[98] 


PCR, ISH 


SCP 


USA 


6, 11, 16, 18 


5/14 


35.7 


Fliederef al. 1998 


[99] 




PCR, ISH 


SCP 


Japan 


6, 16 


1/1 


100 


Haradaef a(. 2000 


[101] 


ISH 


SCP _ 


France 




0/1 


0.0 


Paganin etai. 2009 


[70] ^ 


ISH, PCR 


SCP 


Japan 




0/1 


0.0 


Inamura et al. 2011 


[102] 




SCP 


liSA^ 




0/1 


0.0 


Lagana etai. 2011 


[103] 


ISH 


SCP 


USA 




1/3 


33.3 


Lang et at. 2011 


[74] 



'Method listed first was used as the HPV detection method in meta-analysis database, EM - electron microscopy, IHC - immunohistochemistry, ISH - in situ 
hybridization, PCR - polymerase chain reaction, SCP - squamous cell papilloma 



CI: 0.311-0.542), using tlie fixed effects (FE) model, and 0.495 

(95% CI: 0.316-0.675), using the random effects (RE) mod- 
el (Fig. 1). Table 2 depicts thie meta-analytic results of thiose 
15 studies, stratified by the HPV detection technique. There 
is significant heterogeneity only between the studies (n = 9) 
using ISH, as measured by Cochran's Q statistic, with p = 0.038. 
This heterogeneity becomes of borderline significance in com- 



parison within strata (p = 0.056), and disappears in com- 
parison between strata (random effects model, p = 0.348). 
The percentage of variation (P) is higher (50.9%) for ISH-based 
studies than among PCR studies (7.0%). Using the RE mod- 
el, studies based on ISH give higher point estimates of HPV 
prevalence (0.534, i.e. 53.4%), as compared with PCR stud- 
ies (39.1%). 



Model 



Random 



Study name 


Subgroup 


Method 


Year 




Statistics for each study 








within study 






event 
rate 


lower 
limit 


upper 
limit 


Z-value 


p-value 


total 


Trillo et al. 


SCP 


IHC 


1988 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Kerley et al. 


SCP 


ISH 


1989 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Carey et al. 


SCP 


ISH 


1990 


0.067 


0.009 


0.352 


-2.550 


0.011 


1/15 


Bejui-Thivolet et at. 


SCP 


ISH 


1990b 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Popper et at. 


SCP 


ISH 


1992 


0.929 


0.423 


0.996 


1.748 


0.081 


6/6 


Yousemat et al. 


SCP 


ISH 


1992 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Popper et at. 


SCP 


PCR 


1994 


0.355 


0.209 


0.534 


-1.593 


0.111 


11/31 


Katial et al. 


SCP 


ISH 


1994 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Flieder et al. 


SCP 


PCR 


1998 


0.357 


0.157 


0.624 


-1.054 


0.292 


5/14 


Kawaguchi et al. 


SCP 


PCR 


1999 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Harada et al. 


SCP 


PCR 


2000 


0.833 


0.194 


0.990 


1.039 


0.299 


2/2 


Paganin et al. 


SCP 


ISH 


2009 


0.167 


0.010 


0.806 


-1.039 


0.299 


0/2 


inamura et al. 


SCP 


ISH 


2011 


0.167 


0.010 


0.806 


-1.039 


0.299 


0/2 


Lagana et al. 


SCP 


PCR 


2011 


0.167 


0.010 


0.806 


-1.039 


0.299 


0/2 


Lang et al. 


SCP 


ISH 


2011 


0.333 


0.043 


0.846 


-0.566 


0.571 


1/3 










0.495 


0.316 


0.675 


-0.054 


0.957 


38/89 



Event rate and 95% CI 



Rg. 1. Forest plot of the 15 studies reporting on HPV detection in bronchial SCPs 
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Table 2. Meta-analysis of the 15 studies* stratified by the HPV detection method 



Detection 


No. of Events Sample 


Point estimates of 


Point estimates of 


Homogeneity 


l-squared 


Homogeneity 


method 


studies 




size 




event rates 


event rates 


(Cochran's Q)** 


(|2)** 


(p-value)** 










(fixed effects model) 


(random effects model) 
















Point 


95% CI 


Point 


95% CI 








IHC 


1 


2 


2 


0.833 


0.194-0.990 


0.833 


0.194-0.990 


0.000 


0.000 


1.000 




9 


16 


36 


0.478 


0.268-0.697 


0.534 


0.232-0.812 


16.317 


50.971 


0.038 


PGR 


5 


20 


51 


0.384 


0.259-0.526 


0.391 


0.252-0.549 


4.305 


7.082 


0.366 


Sumnnary 


15 


38 


89 


0.422 


0.311-0.542 


0.434 


0.302-0.575 


22.690 


38.300 


0.065 


Total within (FE) 














20.622 




0.056 


Total between (FE) 














2.069 




0.355 


Total between (RE) 














2.112 




0.348 



'All studies reporting only bronchial carcinomas are omitted; 
"only calculated for fixed effects model 

IHC - immunohistockemistry, ISH - in situ hybridization, PCR - polymerase chain reaction, FE -fixed effects model, RE - random effects model 



When stratified by the geographic origin of the study 
(Table 3), there is significant heterogeneity (p = 0.027) only 
betw/een the studies from Europe, with the percentage of 
variation up to 63.6%. The highest summary effect size 
(63.1%) is derived from studies reported from Other Asia, 
followed by those conducted in USA/Canada (51.7%), and 
Europe (41.1%). In the between-strata comparisons using 
the RE model, the difference in results is not significant 
(p = 0.792), indicating that the heterogeneity between the 
studies from different geographic regions is not statistically 
significant. 

Meta-regression 

In meta-regression testing HPV detection methods as 
study-level covariates (Table 4), none of the methods result- 
ed in summary point estimates that were significantly dif- 
ferent from the reference (IHC), includingthe effect size dif- 
ference (0.052) between ISH and PCR methods {p = 0.626). 
The HPV detection method was not a significant study-lev- 
el covariate (p = 0.150 for regression coefficient pi or effect 



parameter). The same was true when only the studies using 
ISH (n = 9) or PCR {n = 5) were included in this meta-regres- 
sion (p = 0.484). 

Using Other Asia as the reference, the effect size differ- 
ences to the other geographic regions are not statistically sig- 
nificant (p = 0.139 and p = 0.510) (Table 5). In formal meta- 
regression, the geographic origin of the study did not have 
a significant impact on the effect size (p = 0.164). 

Publication bias 

There was little evidence for publication bias among stud- 
ies based on ISH: Begg p = 0.301, Egger's p = 0.019, Duval 
and Tweedie's trim and fill (RE) method identified no miss- 
ing studies. For studies using PCR, there was some evidence 
of publication bias (Begg p = 0.231; Egger's p = 0.161), but 
Duval and Tweedie's trim and fill (FE and RE) method im- 
puted 1 hypothetically missing study, with a marginal effect 
on adjusted point estimates (from 0.391 to 0.373). 

There was no evidence for publication bias among stud- 
ies {n = 3) from Other Asia and Europe (n = 5). For studies 



Table 3. Meta-analysis of the 15 studies* stratified by their geo 


graphic orig 


;in 








Detection No. of Events Sample 
method studies size 


Point estimates of 

event rates 
(fixed effects model) 


Point estimates of 
event rates 
(random effects model) 


Homogeneity 
(Cochran's Q)** 


l-squared 

(|2)** 


Homogeneity 
(p-value)** 








Point 
estimate 


95% CI 


Point 
estimate 


95% CI 








Other Asia 3 


4 


6 


0.631 


0.229-0.908 


0.631 


0.173-0.933 


2.878 


30.510 


0.237 




North America 7 


14 


27 


0.486 


0.297-0.678 


0.517 


0.297-0.730 


6.808 


11.867 


0.339 




38 




0.422 








22.690 ^ 






Total within (FE) 














20.675 




0.055 


Total between (FE) 














2.015 






Ik- 

Total between (RE) 














0.467 




0.792 



'All studies reporting only bronchial carcinomas are omitted 
"Only calculated for fixed effects model 
FE -fixed effects, RE - random effects 
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Table 4. Effect of HPV detection method on the effect size in maximum likelihood meta-regression 



Study-level covariates 


No. of studies 


Effect Size* 


Difference in effect size estimates 


HPV detection method 


(homogeneity p-value)** 


Point estimate 


95% CI 


Difference in 
point estimat^ 


95% CI 


IHC 


1 (p = 1.000) 


0.833 


0.194-0.990 


1.000 






9(d = 0.038) 


0.534 


0.232-0.812 


0.5551 


0.393-0.717 


PCR 


5 (p = 0.366) 


0.391 


0.252-0.549 


0.6072 


0.473-0.741 


Meta-regression for all methods Slope: -0.409 (95% CI: 


-0.968-0.148) (p = 


0.150): Intercept: 2.389 (95% CI: -1.327-6.107) (p = 0.207) 


ISH 


9 (p = 0.038) 


0.534 


0.232-0.812 


1.000 





Meta-regression (ISH/PCR) Slope: -0.386 (95% CI: -1.469-0.697) (p = 0.484): Intercept: 2.231 (95% CI: -5.062-9.525) (p = 0.548) 



'Random effects model 

"Cochran's Q; IHC, immunohistochemistry 

'p = 0.126 
^p = 0.087 
= 0.626 

ISH - in situ hybridization, PCR - polymerase chain reaction. Slope; effect parameter (= regression coefficient /},), Intercept (= coefficient 
Table 5. Effect of geographic origin of the study on the effect size in maximum likelihood meta-regression 



Study-level covariates 


No. of studies 


Effect size* 


Difference in effect size estimates 


Geographic Origin of Study 


(homogeneity p-value)** 


Point estimate 


95% CI 


Difference in 95% CI 
Mint estimates 


Other Asia 


3 (p = 0.237) 


0.631 


0.173-0.933 


1.000 0.321-0.604 


r Europe 


5 (p = 0.027) 


0.411 


0.127-0.770 


0.3091 -0.088-0.707 


North America 


7 (p = 0.339) 


0.517 


0.297-0.730 


0.1482 -0.273-0.569 





'Random effects model 

"Cochran 's Q 
'p = 0.139 
^p = 0.510 

Slope; effect parameter (= regression coefficient pj. Intercept (= coefficient pj 

[n = 7) from North America, there was some evidence for pub- 
lication bias (Begg p = 0.088; Egger's p = 0.0691), and Duval 
and Tweedie's trim and fill (FE and RE) method imputed 
2 hypothetically missing studies, with a marginal effect on 
adjusted point estimates (from 0.517 to 0.441). 

Sensitivity analysis 

Meta-analytic results seemed relatively robust to all {n = 14) 
one-by-one study removals, with little change in the mag- 
nitude and precision of the FE and RE summary point esti- 
mates of the effect size. The single most influential study was 
the one with the highest sample size (n = 31) [98], the removal 
of which would increase the summary effect size from 0.422 
to 0.475 (FE model) and from 0.495 to 0.538 (RE model). 

Discussion 

The role of HPV in etiology of bronchial papillomas has 
received increasing interest since the late 1970s, when the 
first evidence was provided [79-82]. The accumulated lit- 
erature has been reviewed most recently some 10 years ago 
[8, 84], but no formal meta-analysis has been published as 
yet. The formal meta-analysis presented in this communi- 
cation covers all the published literature reporting HPV detec- 
tion in bronchial SCPs {n = 15 studies), and also updates the 



author's own review of 2002 [84]. Importantly, the authors 
did not make any restrictions as to the HPV detection meth- 
ods, although some of the early HPV detection techniques 
(IHC) are obsolete by now. This was just to validate by for- 
mal meta-analysis and meta-regression the frequently pre- 
sented concept that the wide variation in HPV prevalence 
in bronchial SCPs might be explained by different detection 
techniques [8, 84]. The other study-level covariate addressed 
in this meta-analysis and meta-regression is the geograph- 
ic origin of the study, also listed among potential causes of 
variation in HPV prevalence in respiratory tract neoplasia 
[84-86]. 

To evaluate the heterogeneity in meta-analysis is crucial 
because the presence or absence of true heterogeneity (i.e., 
between-study variability) directly affects the statistical mod- 
el that should be used to analyze the database [87, 104-106]. 
The usual way of assessing whether true heterogeneity exists 
has been the Q test, originally introduced by Cochran 
(1954) [107]. Significant p-values in the Q test implicate true 
heterogeneity and warrant the use of a random effects (RE) 
model to test both within- and between-study variability. The 
Q statistic does not indicate the magnitude of true hetero- 
geneity, however, but only its statistical significance [107]. 
On the other hand, the P index, introduced by Higgins etai 
(2003), measures the extent of true heterogeneity, expressed 
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as the percentage of the total between-study variability in 
the effect sizes [108]. One of the major advantages of the 
P index is that the indices obtained from meta-analyses with 
different numbers of studies and different effect metrics are 
directly comparable [104, 108]. 

Given the above considerations, there is little doubt that 
marked heterogeneity exists only between the studies (n = 9) 
using ISH, but not those using PGR [n = 5) methods, as esti- 
mated using the Q test and 1^ index (Table 2). This marked 
heterogeneity among studies comprisingthe bulk (9/15) of 
all included studies justifies the adoption of the RE model 
to analyze the summary statistics for heterogeneity [104-106]. 
Usingthe RE model, the most important conclusion implies 
that there is no true heterogeneity between the studies using 
different HPV detection techniques, as indicated by the non- 
significant homogeneity p-value (p = 0.348) for the between- 
study comparison, also when the EE model is used (p = 0.355). 
In other words, we can revisit the concept raised in some 
recent reviews [8, 84], suggesting that the differences in 
HPV prevalence reported in bronchial papillomas might be 
explained by the different HPV detection techniques. This 
statement is limited, however, by the relatively small num- 
ber of studies reporting on HPV detection in these lesions. 

An alternative view suggests that this variable HPV 
prevalence would be related to the different geographic 
regions of the study origin [8, 84]. To evaluate this concept, 
we performed our meta-analyses stratified by the geographic 
origin of studies (Table 3). Both the Q test and the P index 
demonstrate a statistically significant heterogeneity (p = 
= 0.027, 1^ = 63.6%) only between the studies (n = 5) derived 
from Europe, but not among those from the other two geo- 
graphic regions (Other Asia, North America). However, when 
the RE model was used to calculate the summary statistics, 
the homogeneity p-value was 0.792 for the between-stra- 
ta comparison. This implies that the seemingly wide varia- 
tion (11.8-63.6%) in HPV prevalence from different geographic 
regions is not significant according to the strict meta-ana- 
lytical criteria. Noteworthy, however, is that the highest sum- 
mary effect size estimates (63.1% for Other Asia) are based 
on a small number (n = 3) of studies only, while the lower 
estimates (41.1-51.7%) are derived from at least 5 to 7 stud- 
ies (Table 3). In this respect, we should also keep in mind one 
of the shortcomings of the Q statistic, i.e., it has a poor pow- 
er to detect true heterogeneity in a meta-analysis including 
a small number of studies, but excessive power to detect even 
insignificant variability when a large number of studies is 
included [87, 104-108]. 

We also performed m eta- regress ion to formally compare 
these differences in summary effect sizes [92]. In meta-regres- 
sion with the HPV detection method as the covariate, the 
regression coefficient for the effect parameter (pi, or slope) 
was not statistically significant (p = 0.150). The same is true 
when the geographic origin of the study (p = 0.164) was test- 
ed for its impact as the study-level covariate. These data for- 
mally confirm that HPV detection method and geographic 
origin of the study are not significant study-level covariates 
accounting for the heterogeneity of HPV prevalence in 
bronchial papilloma studies. The results of meta-regression 
did not change if only the studies based on ISH and PGR were 
included in the analysis. Thus, despite the seemingly higher 



summary HPV prevalence derived from ISH studies (53.4%) 
as compared with the PGR studies (39.1%), this difference 
does not reach statistical significance (p = 0.484). 

Some evidence for publication bias was detected, but this 
usually had an insignificant effect on the adjusted point esti- 
mates in stratified meta-analysis. There was no evidence for 
publication bias among ISH-based studies, and the slight pub- 
lication bias revealed among the PGR studies resulted only 
in a negligible drop of the point estimates (from 39.1% to 
37.3%). As to the geographic areas, some publication bias 
was evidenced only among the studies from North Ameri- 
ca Thus, for this region, Duval and Tweedie's trim and fill (both 
EE and RE) method imputed 2 hypothetically missing stud- 
ies, with a marginal impact on summary point estimates 
(51.7% to 44.1%). Thus, it remains to be seen whether the 
future studies from this region will report somewhat lower 
HPV prevalence than the summary effect size (51.7%) 
derived from the currently available studies. 

In sensitivity analysis, based on one-by-one removal of all 
15 studies, all meta-analytic results seemed relatively robust 
to all removals, with no major change in the magnitude or 
precision of eitherthe EE or RE summary point estimates of 
the effect size (Fig. 1). In a meta-analysis based on such a small 
number of studies, and relatively small number of analyzed 
cases (n = 89), however, it is expected that the most influ- 
ential studies are those based on the largest number of cas- 
es, and/or those with very low or very high event rates. In this 
analysis, such a single study was that by Popper et al. [98], 
reporting HPV prevalence of 35.5% (11/31) in a series of 31 
bronchial SGPs. The relative (fixed) weight of this single study 
is 42.9, and even if reporting an effect size (0.355) far below 
the summary effect size (0.495) derived from all studies, 
removal of this study from the meta-analysis would increase 
the summary effect size (random) only from 0.495 to 0.538, 
i.e., by 4.3%. To be truly influential, any such hypothetical study 
should include over a hundred cases and demonstrate HPV 
prevalence markedly deviating (down- or upwards) from the 
summary effect size derived from the existing 15 studies. 

Taken together, these meta-analytical results based on 
all published literature (n = 15 eligible studies) on HPV detec- 
tion in bronchial squamous cell papillomas imply that HPV 
prevalence varies accordingto i) HPV detection method, and 
ii) geographic origin of the study In stratified meta-analysis 
and meta-regression, however, this variability is not signif- 
icantly associated with either of these two study-level 
covariates. Because not formally confirmed by the meta- 
regression, it seems premature to conclude that bronchial 
papillomas in different geographic regions have a different 
etiology, as hypothesized in some recent studies for bronchial 
cancer [109]. Large multi-center studies based on larger series, 
as well as prospective cohort studies, are needed to better 
elucidate the impact of HPV in pathogenesis of bronchial 
papillomas (papillomatosis), which are a not infrequent 
accompaniment of their (synchronous or metachronous) 
malignant counterparts, i.e., bronchial squamous cell car- 
cinoma [2, 8, 11, 27, 39, 40, 50, 51, 57, 58, 66, 71, 84]. 
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